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ABS TRACT 


Seasonal longshore flow patterns are examined at four points along 
the Oregon-Northern California coast. Summer and winter activity is 
examined as far seaward as 25 nautical miles and as deep as 200 meters. 
Long-term mean hydrographic data are used to determine geostrophic ve- 
Mecities. A mearshore baroclinic southward flow (~20 cm/sec) is ob- 
served at each of the points during the summer. Winter currents are 
Gemenrally very small) (<10 cm/sec) and largely barotropic in nature. 
Seasonal volume transports are presented; corrected velocity profiles 
are also presented based on data from moored current meters. Quali- 


tative explanations of the observed phenomena are considered. 
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I. INTRODUCTION 


Pee eoaCRIPTION OF AREA OF STUDY 

The area of the study, shown in Figure 1, is the one studied by 
Bourke (1971). More specifically it is bounded by 41° and 46° North 
latitude. It extends from the coast of Oregon out to approximately 
25 nautical miles. Subsurface activity is examined down to a depth 
of 200 meters. 

The California Current flows through the study area and is the 
main influence on water movement. Significant variations from the 
California Current pattern may develop locally. These deviations 
may be attributed to winds, tides, the shallow topography or a com- 
bination of these elements. 

The topography is characteristic of the West Coast of the United 
States in that the shell? is very narrow and steep. Steepness general- 
ly increases southward as seen in Figure 2. Immediately north of 
Zone A is Astoria Canyon just off the mouth of the Columbia River. 

At the southern end of Zone B are located Heceta Bank and Stonewall 
Bank, features that shoal to depths less than 60 meters and may cause 
some deflection of the longshore current. The protrusion just south 
of Heceta Bank is Cape Blanco which divides Zone C from Zone D. 

A mixture of the Oyashio and the Kuroshio current systems supplies 
Giemealitemtla current, The current is broad, slow, and shallow 
(~500 m). It flows southward in a band roughly 300 miles wide at an 
average speed of 10 cm/sec (Dodimead, Favorite, and Hirano, 1963). A 
seasonal variation of the California Current is the Davidson Current. 


It develops off the Washington-Oregon coast in September, becomes wel] 
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Figure 1. The coastal boundary of the study area. 
(from Bourke, 1971) 
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established by January and disappears by May. It develops inshore of 
tie California Current and may reach 50 miles in breadth and attain 
speeds of 20 to 40 cm/sec (Schwartzlose, 1963). 

The Oregon coast is located approximately in the center of the 
zone of prevailing westerlies with local winds varying from northwest 
to southeast throughout most of the year. The barrier presented by 
the mountains of the Coast Range influences the general wind pattern, 
deflecting the winds so that they tend to align with the trend of the 
coast (Cooper, 1958). In regions where the mountains are low, deflec- 


mone is Minimal and normal wind conditions exist. 


B. ORGANIZATION OF DATA 
The National Oceanographic Data Center (NODC) has on file all hy- 
drographic data taken in this region prior to 1966. Data from 1966 
to 1971 are those observations made by the Department of Oceanography 
of Oregon State University (OSU). The NODC data are filed by 10° Mars- 
den square numbers; number 157 encompasses the area of ane study. 
Bourke employed a computer program which excluded all data taken 
more than 25 miles from shore. The data were then grouped by month so 
that summer and 2 seasons could be determined. The summer season 
included all data sampled in July and August while the winter season 
included data from December through March. These data were found to be 
grouped about hydrographic lines frequently sampled by OSU. These lines 
originate off Astoria (AH line), Tillamook Head (TH line), Depoe Bay 
(DB line), Newport (NH line), Coos Bay (CH line), and off Brookings 
(BH line). Along hydrographic lines observations are made at regular 


intervals at 5, 15, 25, we nautrecal miles ortshore, Occasionally 
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additional stations are taken at three and ten miles. As a result, 
within the study area, most of the data were clustered about hydro- 
graphic lines oriented perpendicular to the coast and were further 
Gueistered at positions 5, 15, and 25 miles from the coast. 

With these points in mind Bourke established a grid system where- 
in average values of temperature and salinity could be computed for a 
specified area within the grid network. Average values were computed 
Meee agepths of 0, 10, 20, 30, 50, 75, 100, 150, and 200 meters. 

The study area was divided into four zones as shown in Figure 3. 
Factors affecting the temperature and salinity were assumed to be uni- 
form or consistent within each zone. The zone boundaries were con- 
structed after having considered the shape of the coastline, location 
of coastal rivers, bottom topography, and distribution of data. Zone 
Peextends from 46° 16' N to 45° 38' N, Zone B from 45° 38' N to 44° 19' 
N, Zone C from 44° 19' N to 42°47'N, and Zone D from 42° 47' N to 41° 
41' N. Each zone encompasses at least one hydrographic line. 

Although the average values apply to a given area, they were as- 
signed to a specific point for calculation purposes. These artificial 
stations provided a network for the computation of geostrophic veloc- 


ie and transport within the study area. 


Cam EREVIOUS STUDIES 

Bourke analyzed the data used in this thesis as regards temperature 
and salinity structure. In general, he found that the surface waters 
were warmer in summer than in winter. Summer temperatures were 2 to 5°C 
warmer than winter temperatures with the largest increases occuring off 


Northern Oregon. Surface temperatures increased with distance offshore 
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Figure 





in both seasons. The offshore temperature gradient was 4 times larger 
in summer than winter. In the longshore direction surface temperature 
increased toward the north in summer while the reverse was true in win- 
ter. The temperature gradient was nearly the same in both seasons. 

Surface waters were 0.5 to 1.0 ppt more saline in summer than in 
winter. Surface salinities decreased with distance from shore in sum- 
mer. In winter surface salinities increased with distance from shore. 
Salinity gradients in the x (offshore) direction were identical in sum- 
mer and winter. In both seasons surface salinities increased toward 
the south. The longshore gradient, however, was twice as large in 
Summer as in winter. 

Variability in temperature, salinity, and density with distance 
from shore was significant in summer but relatively unimportant in 
winter. A strong thermocline was present in summer from the surface to 
30 meters. Winter surface waters were isothermal to 50 meters. A 
strong halocline was observed in both seasons extending to./5 meters 
in the summer and 100 meters in winter. In both seasons at 200 meters 
the temperature and salinity were essentially constant throughout the 
study area. At 200 meters summer temperatures and salinities Wares 
cooler and more saline than winter values due to the strong coastal 
upwelling occuring during the summer months. Summer values were 
peewee? Grand) 33-98 = 0.01 ppt; im winter these were 7.7 £ 0.1°C 
andes oe0 te 0. > ppt. 

Uiemitontat layer (permanent pycnocline), defined by cole 
(1964) to be bounded by the 25.5 and 26.0 o surfaces, rises sharply 


shoreward in summer, but is essentially flat in winter. Its depth in 
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winter is about 100 meters. During the summer it slopes from a depth 

of about 50 meters at 25 miles offshore to a depth of about 20 meters 

at 3 miles offshore. In Zone D it breaks the surface at about 7 miles 
offshore. 

Pavlova (1966) performed a study of the seasonal variations in the 
California Current and Countercurrent in latitudes 20° to 40°N. Using 
hydrographic data averaged over 35 years she computed geostrophic ve- 
locities and transports relative to 1500 decibars. She also mapped 
the dynamic topography of the California Current. 

In developing a characteristic profile of the California Current 
and Countercurrent, she chose a cross section oriented perpendicular 
to the coast around 32° north latitude. A northward countercurrent 
is observed at depth throughout the year. Except during the spring, 

a surface countercyrrent is also observed throughout the year. Max- 
imum countercurrent velocities (8 to 9 cm/sec) are reached in the sum- 
‘mer and autumn at a depth that varies in the course of the year: 

From 200 to 300 meters at the beginning of the summer and in the au- 
tumn, to 100 meters at midsummer. On an annual average, however, 

the countercurrent at depth is practically a single and uniform flow 
moms 25 N to 35 WN and farther. 

The ‘strength of the countercurrent is related to the California 
Current. Maximum velocities for both currents are reached in summer 
and in autumn when the northerly geostrophic wind is strongest and up- 
welling is most intense along the coast. 

Pavlova also notes that in the upper layers there is an upward 


and in the lower layers a downward isopycnic slope from the ocean 
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towards the coast, so that at intermediate depths (100 to 250 meters) 
pmzeone OL Givergence of the isopycnals forms in relation to the posi- 
tion of the profile and the time of year. Her data indicate that the 
countercurrent reaches maximum development in this zone of divergence, 


which in August is centered on the 26.0 isopycnic surface. 


I, 





Ii. THEORY 


A. INTRODUCTION 

It is assumed that the long-term mean current is composed of two 
parts: a baroclinic component arising from the distribution of the in- 
ternal mass field and a barotropic component arising from the slope 
of the sea surface. Geostrophy is assumed in order to determine the 
vertical velocity profile and this is transformed to an absolute pro- 


file by the use of data from moored current meters. 


beeeGtOoTROPHIC VELOCITY 

The theoretical development is presented in the left-handed co- 
ordinate system where the z axis is positive downward and the y axis 
is positive northward. The components of velocity, u and v, are di- 
rected along the x and y axes, respectively. Pressure is denoted by 
Dp, density by p, and the eddy viscosity coefficient by A. Other de- 
Pimitions to be used are specific volume, a, defined as the recipro- 
cal of density and the Coriolis parameter, f, defined as 20sing, where 
@ is latitude and 0 is the angular velocity of the earth. 

Since north-south (meridional) velocities are being examined, 
the appropriate equation of motion describes the east-west (zonal) 


aeceleération: 
do BU Gen ee oe 
dt 0 Ox x 


Higemtrretional term is chosen to have the form: 


lsh Ou 
ath, 5 Ox ee ae +Sca So) 1 (1) 


iipemecdiation states that the frictional forces present can be re- 
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presented as shearing stresses. Furthermore it is assumed that the 
lateral shearing stresses are minor in comparison with the vertical 


shearing stresses. The equation of motion can then be written as 


Qu. 
G 


| 


= ee eau 
oe 


Q. 
CT 


where the V notation has been dropped on the vertical eddy viscosity 
coefficient. 

This equation is simplified by two assumptions based on the fol- 
lowing considerations: 

(a) The data consist of long-term monthly averages of salinity 

and temperature. These values are further averaged over a 
seasonal period of two months or more, so that a smoothed 
seasonal description of the study area may be presented. 

{(b) Frictional effects, although present, are confined to rela- 

tively thin boundary layers (of the order of 20 meters) at 
Ene sea surface and along tne bottom of the shelit and slope 
regimes. ‘These effects exert an insignificant influence on 
Giese costLopiicmilterior.olebile atea. 

In support of (b) the following calculations are presented for 
zone B where sufficient current measurements were available. The 
average wind speed for July and August is 6.0 m/sec blowing roughly 
from the north most of these two months (Fleet Numerical Weather Cen- 
ter, Monterey, California). Based on a wind speed of 6.0 m/sec, Neu- 
man's (1966) value for the "effective'' eddy viscosity coefficient, A, 
is 160 gm/cm-sec. The following current measurements were collected 
by OSU seven nautical miles offshore of Depoe Bay, Oregon, in 100 


metersson water, from 23 June to 10 July 1968: 


My 





Depth of sensor u component* 


2 __(em/sec) 


z= GAS 
50 Zoi 
ue 4.8 


* Velocities are averages based on common record 
Hene@tiie.. (adapted from Pillsbury, et al-., 1970) 


ifaw 1s assumed constant and the above data are used to evaluate 


” (em-sec) 7. 


d° u/d2" numerically, one obtains a value of 2.4 x 107 
The frictional wind stress term is, therefore, of the order of 107° 
males the Coriolis term, fv, is of the order of 10 “. 

By reason of (a) above, du/dt = 0. By reason of (b), 


a A O° u/d2" ~ 0, so that the equation for the x component of motion 


may be expressed as 


meee 
ey: 4 ee 


It is recognized that the time averaging process which eliminates 
the time dependence in equation (2) also reinforces the stress terms 
of equation (1). This result is deemed negligible and equation (3) 
is used as the governing equation for the distribution of the meri- 
dional velocity field. 

For small pressure differences, Ap, and small distances, Az and 


Bxeore may express the vertical and horizontal pressure gradients as 


it 


Ap se aZ 


and Ap (Gay ooo Axe 


it 


Eliminating Ap yields 
g Az = (dp/dx) Ax. 


Rewriting the equation leaves 
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4z _ | op 
6 Ax Ox, 


In the limit that Ax and Az approach zero, while the pressure remains 


constant, one obtains 


dz. _ dp 
BGO ze ox. (4) 


The pressure gradient in equation (3) can now be replaced by the zonal 
slope of an isobaric surface. The concept of a slope, however, implies 
the existence of an absolute reference surface (i.e. a level surface), 
a surface frequently unknown in oceanographic work. Thus the idea of 

a relative slope is introduced. Given the vertical distances, h, and 
hy» between isobaric surfaces at stations A and B, respectively, the 
slope of surface p, relative to a deeper surface ps is given by 


Cy ae 


( 
oh L (5) 


maere | is the distance between stations and station A is east of sta- 


tion B as illustrated in Figure 4. 


B SEA SURFACE A 





V 
Figure 4. Determination of the slope of the p, 
surface relative to the py surface. 
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The evaluation of hy and h, is accomplished utilizing the concept 
of geopotential as described by Sverdrup, Fleming and Johnson (1942). 
The result of their discussion states that a geometrical distance, z, 


is proportional to the pressure integral of the specific volume of the 


water column, 
10 
Lae Gadi 
2 fa ap 


Pmeeemvenient, but arbitrary, choice of units requires that g be spe- 
cified in m/sec”, p in decibars (10° dynes/cm*?) and a in cn’ /gm. The 
depth, z, is then specified in meters. Applying this technique to 


equation (5) yields 


dz 


ae 


10 [ ,\f? Pa ql 
—— On| = Crd 
Gop al % a> Ue 4P,| ¢) 
231 Pi 
Substituting equations (3) and (4) into equation (6) gives the rela- 
tive current perpendicular to the plane between stations A and B, 
_ 10 = Pe Pe ] 
Megs Vou > ey: Cl 7 dp), (j) @ dp), - 7) 
Px P1 

The integral over pressure is the geopotential distance, D, expressed 


imedynamie meters. According to Sverdrup, et al., D is defined as the 


sum of Do and AD, where 


Do 


It 
= 
0 
i 
o 
ee 
—— 
R 
Lo 
wi 
= 
ae) 
ay 
ae 


tf 


Pe 
i Omdp- 
Pi 


and AD 


The standard geopotential distance, Dg, results from integrating over 
a "standard" water column of salinity equal to 35 ppt and temperature 


equal to O°C. The geopotential anomaly, AD, results from integrating 
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over a water column whose salinity and temperature structure are de- 
fined solely in terms of deviations from the structure of the "standard" 
water column. The symbol, 6, denotes the specific volume anomaly. Re- 
writing equation (7) gives 


Pe Pe 
iow =a J Sar, - J dp), |. (8) 
a 1 


Since the ocean is very nearly in hydrostatic equilibrium and 
since the decibar was chosen as the dimension of pressure, dp may be 
replaced by dz, where z is measured in meters. As a working approach, 
a level sea surface is assumed so that v, vanishes. With this assump- 
UO, Vo 1S necessarily a baroclinic velocity only. It is caused 
wentetiy by the distribution of the internal mass field. There is no 
barotropic component since there is no slope of the sea surface. The 
integration is performed from the sea surface to the depth of interest, 


therefore 
Z Z 
v(z) = Vo = ~ (J 6 dz), ‘i (J e dz) |, (9) 
O O 


With this formula v(z) is rendered in meters per second. By speci- 
fying L in kilometers and omitting the factor of 10 in the numerator, 
v(z) may be calculated in centimeters per second. An observer, look- 
menaeethe plane of the stations with station A to the right of sta- 
tion B, will note that a positive current flows away from him while a 
negative current flows toward him. 

By means of a fairly uncomplicated, but lengthy, computer program, 
the specific volume anomaly may be determined by the appropriate poly- 


nomial (Knudsen, 1901). Numerical integration to successive depths 


a 





generates the baroclinic velocity profile. However, the profile is 
still relative to the sea surface which has been assumed level. As 

it turns out this is usually a poor reference level and the calculated 
shear must somehow be fixed at a point to provide a more realistic dis- 
tribution of velocity. Traditionally this has been done by assuming 
that horizontal motion ceases at some great depth, usually around 

1500 decibars (Defant, 1961). ‘The profile is then adjusted to show a 
zero current at 1500 db. In shallower areas a more positive approach 
is to measure the current at some depth (or depths) and then adjust 
the baroclinic profile to make the best fit to the data. In effect, 
this amounts to adding a barotropic component to the velocity profile. 
This latter technique is equivalent to a physical determination of the 
reference level. Both techniques allow the computation of the sea 
surface slope once the shear has been altered. In the computer pro- 
Seams the user Specifies a reference level, and the vertical velocity 
Prot?le is adjusted so that the current is zero at this depth. Know- 
ledge of the correct reference level is imperative if the profiles 


are to provide useful information. 


C. VOLUME TRANSPORT 

When describing oceanic circulation, a quantity of interest is 
the volume transport between two stations. Consider a plane section 
bounded on either side by stations A and B, at the top by the sea sur- 
face and at the bottom by some isobath,,d. The volume transport is a 
measure of the amount of water passing through this plane section per 
unit time. Given the velocity as a function of depth, the volume 


transport is 


Ze 





d 
T= L il VedzZ. (10) 
O 
Using the computer program mentioned previously, it is a simple matter 
to insert the adjusted velocity values into this equation and perform 


the indicated integration. Equation (10) is preferred to the often 


quoted equation, 


10 
ES J @D, SPD Dede (11) 


This equation is developéd by making the following substitution for v 
in equation (10): 


“10 3(AD) _ 10 AD, - AD 


A 
££ Ox f£ 7 


Vv = 


where Dy and D,, are the geopotential anomalies at stations A and B, 
respectively, in the configuration described earlier. Application of 
equation (11) requires that the velocity must vanish at some depth. 

No such requirement is imposed on equation (10). Using the adjusted 
profile, transport in water layers can be calculated and then summed 
HomPive tLOtal volume transport. To a sufficient degree of accuracy, the 


mass transport may be determined by multiplying the volume transport 


by an average density value. 


D. THERMAL WIND RELATIONSHIP 
The thermal wind equation, sometimes written as 


Oh aes 
57 fF Es (12) 


will be referred to in the discussion of procedure. The simplifying 
conditions needed for the development of this expression are given by 


Pillsbury (1972) in a very straight-forward derivation. Using the 
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same coordinate system as before and with only the pressure and Cor- 


iolis forces acting in the east-west direction, the equation of motion 


becomes 
echy 
po fv <a 
Using the hydrostatic equation, 
. 9P 
a8 OZ 


the pressure term is eliminated by cross differentiation to leave 


OM ee Cpe Y Op 
dz fp dx podz. (13) 


In the study area, v is of the order of 10 cm/sec making 


—5 


—= /~~10 
Ras ae) 


aqne slope of the isopycnal, 


we eS . . oo /ox 
OX’ D op /dz 


msenot often less than 10 ~° in the study area (Pillsbury, 1972). 


Therefore the ratio, 


£ OP | v oo 
fp ox p OZ 


is of the order of 10 *%. Thus the second term on the right-hand side 


of (13) may be neglected, and the vertical shear may be written as 


(on) Payee. OP) 
OZ HOM Okes (14) 


Now the static stability is given by 


— 
| 
Ole 
a 
N 


and the product of the isopycnal slope and the static stability is 


24 





given by 


tr 

¢) 

i 

I 
Dlr 
a 


Substituting this expression in (14) leaves equation (12), 


ON eae 
AE £ E Ss. 
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vr EROCEDURE 


A. INTRODUCTION 

The original objective of this study was to compute relative sea- 
sonal velocity profiles in the four zones and then adjust these pro- 
files to agree with observed current measurements. It was thought that 
this procedure would develop a reliable long-term mean picture of the 
water circulation over the Oregon and Northern California shelf-slope 
area. With the corrected velocity profiles, meaningful longshore trans- 
port computations could be made. 

These plans had to be modified due to the lack of sufficient cur- 
rent observations. The best measurements available were collected by 
Oregon State University's Department of Oceanography. These were 
taken in coastal waters off Depoe Bay and Newport, Oregon, in 1966, 1967, 
heos, and 1969. They were used to adjust the velocity profiles of 
Zone B. Other measurements were taken by the University of Washington's 
Department of Oceanography. It was hoped that the latter could be 
Meedeii the profiles of Zone A since they were collected near the mouth 
of the Columbia River. This was not the case, however, since the mea- 
surements that were collected during the appropriate months, were col- 
lected at depths deeper than the available hydrographic data. No 
Current data were located for Zones C and D. In addition scarcity 
of winter hydrographic data in Zone D precluded the formation of sta- 
tistically significant average values for that season. 

Consequently the results presented are not absolute, but do illus- 
Erate a slalrly accurate description of the variations in relative water 


movements from summer to winter and alongshore from north to south. 
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Batis THO 

A computer program was obtained to compute the geostrophic veloci- 
ties via equation (9) as developed in section II. The use of geostro- 
phy is based on the assumption that the frictional boundary layers at 
the bottom and at the sea surface are small in comparison with the depth 
of most of the area. The wind, which is an important consideration in 
this region, is assumed to control only the barotropic component of 
the flow due to its influence on the sea surface slope. 

The problem of wind stress has been reasonably refuted by the nu- 
Menrcal cxample cited in the theoretical development, but some assump- 
tions were involved there that need to be examined. The 1968 current 
data used in the determination of 3°u/dz° were collected in a period 
that did not fall entirely within the months of July and August, the 
period used for geostrophic calculations. It was assumed that during 
these months the currents did not vary significantly from the observed 
Gierents during the time 23 June to 10 July 1968. Also specific wind 
Statistics were not available for the station at the time the current 
measurements were being recorded. Consequently, the wind speed at 
the time of observation was assumed to be the long-term mean wind 
speed as tabulated by FNWC for the months of July and August (July 
Meine ele woeki.s Aucust mean, 11.0 kt.). 

The wind speed was the determining factor in choosing a value for 
the eddy coefficient, A. In order to increase the magnitude of the 
frictional term to the point where it would be significant with re- 
spect to the Coriolis term, it is estimated that the mean wind speed 


from 23 June to 10.July 1968 would have to have been about 20 knots. 





The applicable eddy coefficient would then have a value of about 580 
em/cm-sec instead of the value 160 gm/cm-sec. The frictional term 
would therefore amount to about 10 percent of the Coriolis term and 
probably could not be neglected in treating the equation of motion. 

The assumption that mean wind conditions prevailed at the time 
of data collection might be open to question. Nevertheless, one wishes 
to have confidence in long-term mean data when applying it to any 
one particular time segment. If the number of observations of wind 
speed is large and these observations are concentrated about the mean, 
then the assumption of mean conditions is probably adequate. This 
is the case in the study area during the summer. 

The irregular topography of the area will complicate the geostro- 
phic flow pattern but will not affect the condition of geostrophy 
meGese €xcept im the very shallow regimes. Dynamically the topo- 
graphic anomalies do not appear to influence the hydrographic data. 
For a given depth, temperature and salinity do not vary significantly 
in the longshore direction within the zone boundaries. This fact is 
one of the points considered by Bourke when he first established the 


grid system of this area. 





Ver RESULTS 


A. VELOCITY FIELDS IN SUMMER 

The most pronounced results are presented in Figures 5 through 32. 
These profiles represent geostrophic meridional velocities that were 
computed relative to the deepest common data point between stations. 
ihieyed@ mot represent accurate absolute values, but they are reliable 
iiecetms OL the distribution of vertical shear. Examined in this 
light, the profiles illustrate some interesting patterns. 

The diagram for Zone A during the summer shows a surface jet start- 
ing at 13 to 15 nautical miles offshore and extending to the seaward 
finde or the data (25 n. mi.). This jet is the result of the dis- 
charge from the Columbia River. Using theory developed for surface 
discharges from ocean outfalls, Bourke has shown that the temperature- 
salinity characteristics of myencoré of the Columbia River plume remain 
intact to a distance of 10 to 15 miles from the mouth of the river. 
The jet in Zone A profile, therefore, represents a southward flowing 
mass OL warm and relatively fresh water that has been discharged 
from the jetties at the mouth of the Columbia. The averaged hydro- 
graphic data support this analysis since cold saline oceanic water is 
found on both sides of the jet which is characterized by its warn, 
relatively fresh flow. 

The Columbia River discharge exerts a minimal influence on Zones 
B and C. In these zones the distinctive temperature-salinity of the 
core of the Columbia River plume is not present and oceanic water pre- 
dominates. The Zone B diagram also shows a surface jet located in the 


interior of the zone. Other barotropic data could shift the location 
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Figure 5. Average summer meridional velocity 
field for Zone A in cm/sec. 
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Figure 6. Average summer meridional velocity 
field for Zone B in cm/sec. 
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Figure 7. Average summer meridional velocity 
field for Zone C in cm/sec. 
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Figure 8. Average summer meridional velocity 
field for Zone D in cm/sec. 
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of the jet east or west, but 10 nautical miles is a reasonable posi- 
tion. This jet flow is a characteristic predicted by a numerical model 
of this area developed by O'Brien and Hurlburt (1972). They used 
a model consisting of two homogenous layers coupled together by shear- 
ing stresses but not by mixing. The system is driven by a meridional 
wind stress. Their results predict a southerly current centered 10 
kilometers offshore, in reasonable agreement with the value of 10 
nautical miles (20 km.) predicted here. 

The Zone C profile is based on less data than Zones A and B. Of 
primary concern here is the absence of sufficient hydrographic data 
to construct long-term mean values at a position less than 5 miles 
off the coast. As a result any isotachs drawn shoreward of 10 miles 
are speculation. Nevertheless, the profile indicates that the jet, 
if it is the same one from Zone B, has moved seaward. It is note- 
worthy that here, as in the other zones, the jet remains positioned 
above the steeper portion of the permanent pycnocline as described by 
Bourke. 

What might be happening in Zone C is that the northerly wind 
along the coast is being constricted and accelerated as it passes Cape 
Blanco. The physical result is an area of uniformly high winds and 
decreased pressure. Consequently one observes more intense upwelling 
in this area than in Zone B. 

Zone D, located just south of Cape Blanco, appears to exhibit 
the effect dramatically. Upwelling has intensified as shown by the 
presence of cold saline water farther out to sea and by the permanent 


pycnocline breaking the surface at about 7 miles offshore as mentioned 
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earlier. 

As discussed beforehand these profiles were computed relative 
to the deepest common data point between Stations. One cannot expect 
this method to be very accurate in an absolute sense, but it is more 
realistic than assuming surface currents are zero. As a check on 
the procedure thus far, Figures 9 and 10 are presented. These are pro- 
files for Zone B that have been corrected using current measurements 
jmomemoored Current meters. The data, collected and processed by OSU, 
are listed in Tables I and ITI. 

The measurements from 1966 indicate a vertical shear of 4.7 x 107° 
sec ~~ between the depths of 20 and 60 meters at 5 nautical miles 
offshore. Using mean hydrographic values in the thermal wind rela- 
tion yields a vertical shear of about 4.2 x 10 ~° sec ~* between sta- 
Bions b-J> sand B-l> Over an interval of depth centered at 40 meters. 
The 1969 observations show a vertical shear of 2.6 x 10~° sec * be- 
tween the depths of 20 and 80 meters at 15 nautical miles offshore. 


Between stations B-15 and B-25, the average hydrographic data in the 


thermal wind relation gives a vertical shear of 1.6 x 107° sec 
over an interval of depth centered at 50 meters. The agreement be- 
tween the measured shear and the predicted geostrophic shear is re- 
assuring and lends a little more credence to the adjusted profiles. 
Consistent with the profiles discussed previously, the adjusted 
profiles (Figures 9 and 10) confirm the existence of the surface jet. 
The presence of a subsurface countercurrent is indicated by the 1966 
profile. These results have also been noted by Mooers, Collins, and 


Soe) eee incin data imeluded that of Table I and a line of 
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Figure 9, Summer meridional velocity field for Zone B 
adjusted with 1966 current observations. 
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Figure 10. Summer meridional velocity field for Zone B 
adjusted with 1969 current observations. 
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japle ies Mean Current Statistics, 1966 


(based on hourly averages for common 
record length of 339 hours) 


Depoe Distance Water Sensor Eastward Northward 
Bay offshore depth depth component** | component** 
Station (n mi) (m) (m) (cm/sec) (cm/sec) 
80 20 —2.2# 11 Sie Sa 1 
80 60 2.8 = 6 Oy aie, 
10 10 140 20 is ie ee, “12.1 +9 
eS L5 200 60 6.0 + 6 fe9ee 13 


* Adapted from Mooers, Collins, and Smith (1972). 


** Velocities are presented as the mean value + the standard 
deviation. 


Table I1.* Mean Current Statistics, 1969 


(based on total record lengths 
acquired in August and September) 


Station Distance Water Sensor Eastward Northward 
(Newport or offshore depth depth companent** component** 
Depoe Bay) (n mi) (m) (m) (cm/sec) (cm/sec) 

NH-3 3 50 20 O.7 + 7 “4.4 4+ 16 
DB-/ 7 100 40 -6.1 + 11 0 | ee sel AS) 
NH-15 15 100 20 “9.9 + 16 @16 103 16 
NH-15 15 100 80 ar On “2.8 + 12 


* Adapted from Huyer, et al. (1971). 


Velocities are presented as the mean value + the standard 
deviation. 
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hydrographic stations taken in late July, 1966, off Newport, Oregon. 
The velocity field presented by Mooers, et al. differs only slightly 
from the profile presented here based on long-term mean hydrographic 
data and the data from Table I. The 1966 profile centers the south- 
ward flowing jet at about 5 nautical miles (10 km.), in better 
agreement with the results of O'Brien and Hurlburt. The northward 
undercurrent apparently does not break the surface as suggested in 
Meeure 9. ~Mocers, et al: indicate that in 1966 the surface flow 

was all southerly at least as far as 60 miles offshore. They report 
that the undercurrent may have come within about 10 meters of the sur- 


Pace that summer. 


Bee VELOCITY FIELDS IN WINTER 

The winter profiles are presented in Figures 11 through 13. 
Prior to assuming a reference level at depth, the relative (baroclin- 
ic) current values were observed to be very small, illustrating the 
importance of the barotropic component during winter. A zonal sea 
surface slope of one millimeter per kilometer, corresponding to a cur- 
rent speed of the order of 10 cm/sec, would be sufficient to reverse 
the flow at any point in any of the zones. 

In spite of such small geostrophic velocities a description of 
a characteristic flow pattern is attempted. It is felt that the long- 
term nature of the data used in this analysis provides a firmer 
Sround Shomesuch atieattempt. 

Nearshore surface waters in Zones A and B are responding to the 
wind stress while the offshore and deeper waters are moving under the 


influence of the large scale oceanic circulation. Current velocities 
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Figure ll. Average winter meridional velocity 
field for Zone A in cm/sec. 
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Figure 12. Average winter meridional velocity 
field for Zone B in cm/sec. 
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Figure 13. Average winter meridional velocity 
field for Zone C in cm/sec. 
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are greatly reduced from the summer values due to the nearly level 
isopycnal surfaces. Runoff and water piled up onshore by the wind 
contribute a northerly barotropic component to the circulation. 

The University of Washington has made a few current observations 
during the winter season at moorings near the Astoria Canyon (Hopkins, 
1971). Unfortunately, these observations were not sufficiently dis- 
tributed spatially to allow an absolute rendition of the velocity 
field. In each of these cases, however, the mean northerly component 
was only slightly larger (in the sense of being more positive) than 
the predicted geostrophic component. This fact illustrates two points: 
that the fresh water runoff and the water piled up onshore by the pre- 
vailing wind govern the barotropic component of the velocity and that 


this barotropic component is not very large. 


C,. VOLUME TRANSPORT 

Volume transports were computed using equation (10). As in the 
velocity computations, the deepest depth of the data common to the 
EvOustations was used as the reference level. Again it 1S cautioned 
that this is not necessarily the correct level of no horizontal motion, 
but it is better than a level sea surface approximation. 

The results of this endeavor are presented in Tables III and IV. 
As expected, the transports are uniformly small, the winter transports 
being roughly an order of magnitude smaller than the summer transports. 
(This observation is probably correct; however, it is not completely 
justified from an inspection of the tables alone. The winter data are 
generally shallower than the summer data, thereby eroding a basis for 


comparison.) Reflecting the structure of the velocity fields, the trans- 
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Table III. Geostrophic Meridional Transports (Summer) 


Stations 


A-03/A-05 
A-05/A-15 
A-15/A-25 
B-03/B-05 
B-05/B-15 
B-15/B-25 
c-05/cC-15 
C-15/c-25 


D-05/D-15 
D-15/D-25 


Volume Transport 


(Sverdrup) 


0.010 
0.067 
0.411 


O0US 
O20) 
Orel g3 


Orth? 
CesT 


0.134 
0.065 


Direction 


Reference Level 
(m) 
50 
15 
200 


75 
150 
200 


75 
200 


100 
200 


Table IV. Geostrophic Meridional Transports (Winter) 


Stations 


A-03/A-05 
A-05/A-15 
A-15/A-25 


B-03/B-05 
B-05/B-15 
B-15/B-25 


C-05/C-15 
C-15/C-25 


Volume Transport 


(Sverdrup) 


0.001 
0.006 
0.198 


0.001 
02002 
6; 023 


0.034 
0.043 
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Direction 


Reference Level 
(m) 
S10 
(es 
200 


50 
75 
is 


200 





port appears to be southward in both summer and winter. While this 
result is probably true for the summer months, it is uncertain that 
this is the case during the winter months. Deeper data and an accu- 
rate determination of a reference level could drastically change the 
winter transports. The small nearshore northward transport observed 
in Zones A and B during the winter shows the effect of the wind in 


moving the surface waters. 
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V. CONCLUSIONS 


The fact that long-term data were used in this discussion allows 
more general statements to be made than could be made otherwise. 
The surface jet, for example, is shown to be a regular phenomenon of 
the summer season. Its position above the steepest portion of the up- 
turned pycnocline and its absence in winter tie it to the occurrence 
of coastal upwelling. The fact that it is a dominant feature of the 
average summer velocity fields indicates that it must be a steady 
state event or at least an average of a series of consistent events. 
Furthermore, the jet is shown to be a characteristic of most of the 
Oregon and Northern California coastline. 

Mooers, et al. have shown that a two-layered baroclinic model 
of the flow over the Oregon shelf and slope is at least partially 
valid. They postulate that Prem baweteee ie and baroclinic radii of 
deformation should provide fundamental length scales. The barotro- 
pic radius of deformation is defined as the distance a shallow water 
wave would travel in an inertial period: R,__. | a (eit)? f-" , 

barotropic 

The baroclinic radius of deformation is defined as the distance an 
internal wave would travel along the interface between a shallow 


upper layer and a deep lower layer of water in the same inertial 


1 
aa z % ood ee 
period: ee (gD Ap/p)* ££". The estimated values are 
R .. = 2000 kilometers and R . .. = 20 kilometers. By 
barotropic baroclinic 


plotting the depth of the pycnocline as a function of offshore dis- 
EincewOnescMi-Vor-arithmic axes, Mooers, et al. confirmed the expo- 
nential character of the "interface'' and determined its radius of 


jenornatrtonetompomabout G0 kilometers over the shelf and 70 kilometers 
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Svenvetnie slopes by comparison then, the water is deemed sufficiently 
stratified to warrant use of a two-layered (baroclinic) model, al- 
though different length scales are in effect over the shelf and the 
slope regimes. 

Without becoming involved in cause and effect arguments, it is 
asserted that the longshore northerly wind, acting on a relatively 
thin surface layer, moves water offshore. This lighter surface 
water moves seaward without significantly disturbing the lower hea- 
vier water and causes a pressure gradient force directed onshore, 
which in turn gives rise to a southward flow. Because the southerly 
flow diminishes seaward, vorticity is decreased (f = constant, 
ov/dx < 0). If potential vorticity is to be conserved (Stommel, 1966), 
then the thickness of the surface layer must decrease in the vicinity 
of the horizontal shear. Thus one observes the surface jet, the 
sloping pycnocline, and upwelling as interrelated phenomena. 

Another feature that appears to be related to the upwelling pro- 
cess is the poleward flowing undercurrent. The presence of this sub- 
surface flow has been reported by Mooers, et al. and by other workers, 
but no theory has been developed to explain it. It is undeniably 
present in the 1966 profile of Zone B (Figure 9) and, given more 
shallow reference levels for geostrophic computations, it is not dif- 
ficult to imagine its presence at the seaward edge of the profiles 
presented in Figures 5 through 8. 

It is reasonable to expect that the inclined isopycnals nearshore 
are capable of producing an offshore pressure gradient force if the 


sea surface does not develop enough slope to counteract it. Under 
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the proper conditions therefore, one may obServe the resulting north- 
ward baroclinic flow located beneath the pycnocline and offshore. In 
the shallow shelf regime, the effect is masked by shearing stresses 
with the bottom and with the nearshore surface jet. This may explain 
why the subsurface flow is found during some summers and not found 
during others. 

These conclusions are compatible with those of Mooers, et al. They 
are in basic agreement with the results of O'Brien and Hurlburt, but 
here there is a notable discrepancy. Upwelling is known to consist of 
a tranSient and a steady state condition. This discussion cannot pro- 
vide an insight into any transient solutions, but it does indicate 
that the nearshore jet and possibly the subsurface flow develop into 
steady state features. O'Brien and Hurlburt, in obtaining the near- 
shore jet, do not permit this result. Indeed their model does not 
have a steady state solution so that a discussion of such fundamental 


da@rrercnces may not be valid. 
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VI. SUMMARY 


The features observed in the study area are attributable to gen- 
eral oceanic circulation and to local wind and climate. The California 
Current is the main influence on water movement over the shelf and 
slope. In general, the flow is southward regardless of the season. 
Summer winds are northwesterly while winter winds are southeasterly. 
These winds govern the sea surface slope and thereby determine the 
barotropic flow during the different seasons. Surface waters re- 
spond to the wind field and, in general, flow parallel to it. Cli- 
mate along the Oregon coast provides plentiful fresh water runoff 
during the winter and considerably less runoff during the summer. 

The influence of the Columbia River discharge that is felt so strongly 
in Zone A during the summer is much less evident during the winter 
due mainly to the shift in wind. 

As a fresult of these characteristics more activity is observed 
in the summer than in the winter. The summer wind pattern stimulates 
the process of upwelling which causes tilting of the isopycnals which 
in turn generates the baroclinic flow. The winter wind pattern piles 
up water onshore but apparently to Be cote extent. Consequently 
there is little barotropic flow and even less baroclinic flow. 

Summer circulation is characterized by the upwelling of cold 
saline water next to the coastline. Typically the surface salinity 
increases and the surface temperature decreases shoreward and south- 
ward indicating more intense upwelling in the lower latitudes of the 
area. Also observed during the summer is a nearshore surface jet 


flowing southward. The surface jet occurs in the region above the 
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steepest portion of the permanent pycnocline where the near surface 
isopycnals have relatively small slope. In Zone A the jet appears 
to be caused by the discharge of the Columbia River, whereas in the 
other three zones it appears to be intimately related to the upwelling 
process. A subsurface poleward flow is observed but is not conclu- 
sively documented as a regular occurrence. This flow develops below 
the pycnocline in a region where the isopycnals are tilted sharply 
under the influence of upwelling. 

Winter flow is characterized by small barotropic and baroclinic 
components. Surface salinities and temperatures increase seaward 
and southward indicating the greater effect of runoff in the higher 
latitudes of the study area. The very small velocities presented 
in Figures 1] through 13 suggest that no distinctive longshore flow 
pattern exists during the winter season, although the surface waters 
in Zones A and B appear to be responding to the wind field. 

In both summer and winter the volume transport is mostly south- 
ward. Winter transports are consistently less than summer transports 


and they are usually less by an order of magnitude. 
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